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The zero-temperature magnetization process of a one-dimensional S = | spin system with next- 
nearest-neighbour coupling a and anisotropy 8 is studied. There is a region in the parameter 
space where the magnetization has a jump at certain external field. The boundary of this region 
is determined. The two cases are distinguished according as a < j or a > j. In the latter case, 
the energy dispersion of the low-lying excitations exhibits a double minimum structure. 
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The magnetization process is one of the interesting 
topics in the study of low-dimensional quantum spin sys- 
tems. In particular, the one-dimensional systems at low- 
temperature attracted a great deal of theoretical and ex- 
perimental interest because of pronounced effect induced 
by quantum fluctuation in those systems. They ex- 
hibit various nonlrivial behaviours, e.g., the middle-field 
cusp singularity the magnetization plateavfitf and the 
second-order spin-flop transition.^ What we concern in 
the present study is the magnetization jump similar to 
the metamagnetic transition in the classical Ising-like 
antiferromagnet, but it is caused by a different mech- 
anism. 

In this paper, we study a spin-i XX Z chain with 
nearest-neighbour (NN) coupling J and next-nearest- 
neighbour (NNN) coupling Ja in a magnetic field h. The 
model is described by the Hamiltonian 
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where Si's represent the S = \ spin operators on the 
chain, and L is the total number of spins. The periodic 
boundary condition Si+L = Si is assumed. Hereafter, 
we set J > and take a unit so that J = 1. The model 
has been studied extensively. In the absence of exter- 
nal fields, it exhibits a variety of phenomena associated 
with cmimfiting interaction such as spontaneous dimer- 
izatioroBCP and incommensurate spin correlations nu 

We shall concentrate on the region < a < 1 < 
S < where an anomalousJaehaviour in the magnetiza- 
tion process was reported:t3> In a region of parameter 
space, namely, Sf(a) < S < 5 c (a), there exists a critical 
field h c at which the magnetization jumps to its max- 
imum value m s = i (see Fig. 1). It should be noted 
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Fig. 1. The zero-temperature magnetization curves for a = 0.25 
and 8 = 0.0, —0.25, —0.5. The thin lines show the finite-size 
results for L = 18 and 20. The magnetization (per site) m is 
normalized by the saturation magnetization m 3 = i . The thick 
lines represent the estimated limiting (L — » oo) curves. For 
8 = —0.25 and —0.5, there exists a critical field h c at which the 
magnetization jumps to the saturation value m s . 



that the model has competing NN and NNN exchange 
couplings in the easy-plane (Xy-plane), while there is 
no competition in the z-axis although couplings are rela- 
tively weak. It is responsible for the existence of the mag- 
netization jump: The system prefers to oder in the XY- 
plane but NNN coupling introduces frustration. There- 
fore, it becomes energetically preferable to order along 
the z-axis for sufficiently large a and \S\. For 6 < 5f, the 
zero-field ground state becomes the completely aligned 
(ferromagnetic) state with m = |. 

The purpose of this paper is to determine the bound- 
ary 5 — S c (a) of— the region where the magnetization 
jump takes place. 13 

First, we discuss the condition that the magnetization 
jump occurs. For the magnetization curve to be smooth, 
eo(m), the lowest energy eigenvalue (per site) for the 
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Fig. 2. (a) Schematic curve of eo(m) in the case when the magne- 
tization jump occurs. There is a region in which £q is negative. 
The line with slope h c is also shown where h c is a critical field, 
(b) Magnetization curve corresponding to eg of (a). 



Fig. 3. (a) The m-dependence of eo and (b) magnetization curves 
of the S = ^ XX Z chain with NNN interaction in the case when 
Sf < S < S c . The dotted line indicates the critical magnetization 
m c (S). 



Hamiltonian H in the fixed m = J2&i/L sector, must 

be monotonically increasing in m and convex for < 
l . 

2' 



m < m s 
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On the other hand, if there exists a region where Cq < 
0, then the magnetization jumps over that region as 
schematically shown in Figs. 2(a) and 2(b). Thus, the 
presence of such a region as Cq < means that the mag- 
netization jump occurs. The quantity £o( m ) can be es- 
timated from 



£q(to) = lim L x A(m), 



(5) 



where 



A(m) = E (m + + E (m - - 2E(m). 

In the case of our model, the problem of determining the 
boundary curve <5 = <5 c (a) is simplified. The results of 
the exact diagonalization suggest that the critical mag- 
netization m c increases as S is increased and reaches m s 
at S = S c . It indicates that 



(6) 



at least for 5 < S c [see Figs. 1, 3(a) and 3(b)]. The bound- 
ary S c is a line where eg| m= i_o changes sign. Therefore, 
we have merely to consider up to two-spin deviations 
from the completely aligned state for our purpose. We 



put A« = A(± - 1); 



A (D 



where Eq ^ denotes the energy of the ferromagnetic state 
and Eq the lowest energy of the states with n-spin de- 
viations from the ferromagnetic state. For S > Sf, the 
ferromagnetic state is not a ground state of the Hamilto- 
nian H n . However, we shall call those n-deviation states 
as n-magnon states although they are not the excitations 
from the ground state. 

We now obtain the energy spectrum of one- and two- 
magnon states. We rewrite the Hamiltonian (2) in term 
of fcrmion creation and annihilation operators with mo- 
mentum k, a~£ and afc, via Jordan- Wigner transforma- 
tion. The Hamiltonian becomes 
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where 



p(0) 

Wo 



= \{l + a)5L 
l.,, = cos 9 — 5 + a(cos2q — S) 
Ti, 2,3,4 = Scos(q 2 - q 4 ) + aS cos (2(q 2 - q 4 )) 
+2acos(2(7i + q 2 - 93) 

and we have used shorthand notations 1,2,3,4 for 
Qi, 12, <?3, Qi- In cq. (8), Eq ^ and Lu q represent the en- 
ergy of the ferromagnetic state |0) and the excitation 
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Fig. 4. Energy spectrum of one-magnon states. For a > j, two 
minima appears at q = ±qo where go = arccos(— ■£-)■ 



spectrum of the one-magnon state with momentum q, re- 
spectively. For a < j, uj q has a single minimum at q = n, 
while, for a > j, the minimum shifts from q = ir and 
two minima appears at q — ±(jo where qo = arccos(— j^) 
(Fig. 4). 

The total wave number Q of two magnons to be con- 
served, thus we write two-magnon state as 



where 



g>0 



Q = qi+q 2 , q= ~(qi - q-i) 



(9) 



Let us solve the eigenvalue equation 
H\Q) = E\Q), 



(10) 



where E is the energy eigenvalue of the two-magnon 
state. Taking the inner product of both side of eq. ( |Io| ) 
with every basis state {0\a qi a q2 , we obtain the following 
equation: 

4 x - 
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Here, 7q i9 represents the energy of the free two-magnons 
with total wave number Q = q\ + q 2 , and forms a con- 
tinuum for each value of Q. 

We are interested in whether there exists a solution 
of eq. (11) which gives A^ 1 ) < 0. Hence, it is sufficient 
to look for the solution satisfies E < 7q, 9 , i.e., the two- 
magnon bound states below the continuum [cf. eq. (0)]. 
Then, the eigenvalue equation reduces to (see, for in- 
stance, ref. n~3) 
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Fig. 5. The results of AW for a = 0.25 and a = 0.3. There is a 
point <5 C where A' 1 ) vanishes, which is a boundary of the region 
where the magnetization jump occurs. 
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It is the sufficient and necessary condition for the bound 
state lying outside of the continuum to exist; i.e., such 
bound states exists if and only if eq. (|l2|) folds. The 
explicit solution of ( p"2| ) is found analytically for several 
values of a, S and Q. It is investigated numerically for 
general case. 

Figure 5 shows the results of for a = 0.25 and 
0.3. It can be seen that there exists a point 5 C at which 
A^ 1 ) vanishes for each value of a and that A^ is neg- 
ative for iL< 6 C which suggests the magnetization jump 
to occur.H In Fig. 6, the result is compared with the 
one obtained from the diagonalization of the finite chain 
Hamiltonian for rather large systems (L = 20-156) at 
a = 0.25, S = —0.25. Note that there is a weak oscil- 
latory behaviour in the finite-size results. For a > j, 
there appears rapid oscillation, thus it becomes difficult 
to estimate a limiting value of A^ from the finite chain 
calculation. 

In Fig. 7, the curve S c (a) is shown in the a-S plane. 
The magnetization jump occurs in the region enclosed by 
the two curves 5 C and <5/li3^ The ground-state phases in 
the absence of external fields are also shownjtZl 1 there 
are three phases, namely, the spin fluid phase which 
is characterized by the gapless excitation and critical 
behaviour, the spontaneously dimerized phase and the 
ferromagnetic phase. From Fig. 7, it can be clearly seen 
that the curve S c has a kink at a = \. It results from the 
fact that the dispersion curve of the one-magnon states 
changes its shape at a = \: For a < 4, it has a sin- 
gle minimum at q — ir, while, for a > j, two minima 
appears. The lowest lying state should consist of bound 
pairs of magnons with q ~ ir for the former case and with 
±q for the latter case. 
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discussion and comments. 
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Fig. 6. Comparison with the result from the numerical diago- 
nalization. Open circles denotes the finite chain results for 
L = 20-156. The result obtained from eq. (jij) is A' 1 ) = 
-1.43 X 10" 5 . 





-0.2 
-0.4 
<<> -0.6 
-0.8 



-1+ 2 

-1.2 



Spin fluid 



T 1 

\ Dimer 



a o 



_L 



Ferro 



0.125 0.25 0.375 0.5 

a 



[1] K. Okunishi, Y. Hieida and Y. Akutsu: Phys. Rev. B 60 

(1999) R6953. 
[2] K. Hida: J. Phys. Soc. Jpn. 63 (1994) 2359. 
[3] K. Okamoto: Solid State Commun. 98 (1995) 245. 
[4] T. Sakai and M. Takahashi: Phys. Rev. B 60 (1999) 7295. 
[5] F. D. M. Haldane: Phys. Rev. B 25 (1982) 4925. 
[6] K. Okamoto and K. Nomura: Phys. Lett. A 169 (1992) 433. 
[7] K. Nomura and K. Okamoto: J. Phys. A 27 (1994) 5773. 
[8] T. Tonegawa and I. Harada: J. Phys. Soc. Jpn. 56 (1987) 

2153. 

[9] S. R. White and I. Affleck: Phys. Rev. B 54 (1996) 9862. 
[10] C. Gerhardt, K.-H. Mutter and H. Kroger: Phys. Rev. B 57 
(1998) 11504. 

[11] In ref. |lo| Gerhardt et al. determined the boundary <5 C by use 
of the numerical diagonalization of finite chain Hamiltonian 
up to L = 50. However, their result suffers from finite size 
effects and somewhat suspicious (compare Fig. 6 in this work 
with FIG. 3 of ref. [u]). 

[12] I. Ono, S. Mikado and T. Oguchi: J. Phys. Soc. Jpn. 30 (1971) 
358. 

[13] The integrand of (13) has, however, a sharp peak which be- 
comes singular as E approaches to the lower edge of the scat- 
tering continuum 7Q, 9 . Thus, for npenerical quadrature we 
use the double exponential formulajiff which is one of the 
most optimal and suitable numerical quadrature formula for 
improper integrals. 

[14] H. Takahashi and M. Mori: Publ. RIMS, Kyoto Univ. 9 (1974) 
721. 

[15] To discuss thermodynamic limit, we must take into account 
the finite size correction: If A' 1 ' behaves as A' 1 ' = A^ + 
ciL -1 + 0(L~ 2 ) and if A^ = 0, then it is the coefficient 
of the correction term ci that determine the sign of e ' in 
thermodynamic limit. We are unable to estimate the finite 
size correction analytically. However, it seems that there is no 
significant difference in <5 C . 
[16] The boundary of the ferromagnetic phase 8f was determined 
from the level crossing of the ferromagnetic state and the low- 
est eigenstate with m = (see ref. |l7| for deta ils). 
S. Hirata and K. Nomura: :orid-mat70909290| . 



[17] 



Fig. 7. The phase diagram in the a-S plane. The magnetization 
jump occurs in the region enclosed by the two curves <5 C (o) and 
8f (•). The zero-field ground-state phase diagram is also shown; 
the boundary of the ferromagnetic phase 8f and the dimer-spin 
fluid phase boundary (dashed line). 



To summarize, we determined the boundary of the 
region in which the magnetization jump occurs in the 
S = | XXZ chain with easy-plane anisotropy and NNN 
interaction in the case of < a < |. It was shown that 
the a-dependence of the boundary curve S c is qualita- 
tively different according as a < j or a > j. For a > j, 
the region that the magnetization jump occurs shrinks 
with increasing a in contrast to the case of a < 4. It 
reflects the fact that the energy dispersion curve of the 
low-lying excitation develops the double minimum struc- 
ture beyond the line a = i. 

One might be interested in how the system is in the 
case when the external field h is close to the critical field 
h c . One of the possible situation is that there are clusters 
of up-spins separated by fluid- like domains. 
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